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A new data analysis technique for the identification of static and dynamic aerodynamic 
stability coefficients from wind tunnel test video data is presented. This new technique was 
applied to video data obtained during a parachute wind tunnel test program conducted in 
support of the Mars Exploration Rover Mission. Total angle-of-attack data obtained from 
video images were used to determine the static pitching moment curve of the parachute. 
During the original wind tunnel test program the static pitching moment curve had been 
determined by forcing the parachute to a specific total angle-of-attack and measuring the 
forces generated. It is shown with the new technique that this parachute, when free to 
rotate, trims at an angle-of-attack two degrees lower than was measured during the forced- 
angle tests. An attempt was also made to extract pitch damping information from the 
video data. Results suggest that the parachute is dynamically unstable at the static trim 
point and tends to become dynamically stable away from the trim point. These trends are 
in agreement with limit-cycle-like behavior observed in the video. However, the chaotic 
motion of the parachute produced results with large uncertainty bands. 


Nomenclature 


C m 

Pitching moment coefficient 

m. a 

Apparent mass, kg 

Cm a 

Pitch moment slope, 

n 

Number of data points 

Cm & 

Pitch damping coeff., — — 

P 

Pixels 

Cm a 

Linearized C m intercept 

q 

Dynamic pressure, N/m 

D 0 

Reference diameter, to 

R 

Parachute length, to 

9 

Gravitational acceleration, ( m/s 1 2 ) 

Re 

Reynolds number 

h 

Vertical height, m 

So 

Reference area, m 2 

*yy 

Moment-of-inertia, kg — to 2 

t 

Time, sec 

K 

Pixel-per-meter constant 

V 

Velocity, m/s 

l 

Streamwise length, m 

x, y, z 

Cartesian coordinates 

M 

Mach number 


Angle-of-attack, sideslip 

M 

Static moment, N — in 

Ot-T 

Total angle-of-attack 

TO 

Parachute mass, kg 

9, ip 

Pitch, yaw angles 
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p 

Density, kg/m 3 

y,z 

Horizontal, vertical pixels 

<t> 

Clock angle 

00 

Freestream 

Subscripts 

1,2,3 

Calibration point indices 

Superscripts 

/ 

Calibration heights, m 

a 

Apex 



b 

Backshell 

Acronyms 


c 

Camera 

DGB 

Disk-gap- band 

cp 

Center-of- pressure 

MER 

Mars Exploration Rovers 

s 

Swivel 

MPF 

Mars Pathfinder 

p 

Pivot 

TDT 

Transonic Dynamics Tunnel 


I. Introduction 

For the Mars Exploration Rover (MER) mission and other planetary missions using axisymmetric parachutes, 
the stability and flight characteristics of the parachutes can have a critical impact on terminal descent events 
prior to landing. The computer algorithms taking data from accelerometer and gyroscopes to control entry 
descent and landing events (e.g., retro rocket firing) can be adversely affected by erratic motions induced by 
parachute motions. It is therefore important to understand parachute aerodynamics sufficiently to model 
this erratic behavior in simulations prior to flight in order to evaluate the robustness of the integrated landing 
systems in the terminal descent (i.e., under parachute) phase of such missions. 

Historically, it has been a very difficult problem to accurately model the behavior of axisymmetric parachutes. 
While the equations of motion for non-rigid parachutes have been developed, including their dynamic stabil- 
ity, 1 ’ 2 extracting the aerodynamic coefficients experimentally has proven difficult. The aerodynamics from 
rigid models have provided general characteristics of parachute behavior and sensitivities to geometry. 3 How- 
ever, these tests could not address the effects of fabric permeability or non-rigidity of the parachute on the 
aerodynamics. 

Recently, Cruz et al. 4 tested scale versions of the Viking and Mars Pathfinder disk-gap-band (DGB) 
parachutes and variants of those parachutes in the NASA Langley Transonic Dynamics Tunnel (TDT). 
Steady state drag measurements were made with the parachute attached to a swivel at the back of a model 
MER capsule backshell mounted on a strut in the tunnel in a subsequent test. This provided drag infor- 
mation, but could not explicitly measure stability characteristics. The parachutes were then held fixed in a 
different test rig and rotated through a sweep of angles-of- attack to extract static pitching moment, normal 
force and axial force data. While a flight-like fabric parachute was tested, the effects of artificially holding 
the parachute at a specific angle-of-attack on the static stability measurements could not be quantified. The 
work presented here attempts to quantify those effects by extracting the pitching moment of one of the 
parachutes reported by Cruz 4 while in the drag test (i.e., free-to-rotate) configuration. 

Video was obtained during the free-to-rotate and static tests. Cameras were placed immediately downstream 
of the parachute looking upstream, and to the side looking across the tunnel. Knowledge of the dimensions of 
the parachute, tunnel and sting apparatus, as well as the downstream camera position, permits the calculation 
of angle-of-attack, a, and angle-of-sideslip, fi, from the downstream video data of the free-to-rotate tests. 
This a, (3 information can then be used to extract aerodynamic characteristics of the parachutes tested by 
the technique described herein. 

To better understand and model the motion of DGB parachutes, the video of the free-to-rotate test was used 
in an attempt to extract quasi-static aerodynamics for comparison with the static data obtained by Cruz. 
The extraction of dynamic stability derivatives was also performed. This data reduction was performed as a 
proof-of-concept exercise to develop the technique. The test setup and instrumentation was not selected or 
designed specifically for this type of parameter identification. It was hoped that static aerodynamics might 
be extracted with some confidence and evidence of pitch damping functionality with total angle-of-attack 
might be identified. 
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II. Test Setup 


Tunnel Conditions 

The test setup used by Cruz 4 and for the present analysis is shown in Figure 1. Tests were conducted in 
the NASA TDT. All tests were conducted in air at ambient pressure. For the data set presented here, the 
freestream Mach number was 0.102 and the dynamic pressure was 736 N/m 2 (15.4 psf). Shown in Figure 
1 are the two test configurations. The first setup (Figure la) allowed the parachute to move freely in the 
tunnel and was intended to yield drag data only. Video of the parachute in this configuration is the source 
of the data set used in the present analysis. The second setup (Figure lb) held the parachute fixed in a test 
rig. The rig, extending from a turnable in the tunnel wall holds a parachute at its swivel point as well as 
its apex. This allows the parachute to be held at a fixed angle-of-attack (adjustable via the turntable) to 
measure pitching moment, normal force and axial force across a range of angles-of-attack. A ring at the apex 
permits the canopy to slide in the axial direction as the parachute geometry changes with tunnel conditions 
and attitude. Video data of one of these tests was used for a calibration of the data extracted from the drag 
test video. 

The parachute used in this work is a scale version of the Mars Pathfinder (MPF) parachute, although material 
differs from the MPF flight parachute. The test parachute used in this data analysis is Fill Nylon, which 
is less permeable than the fabric used on the MPF flight parachute. The dimensions and flow conditions of 
the test are listed in Table 1. 

The coordinate system used in the test setup is as follows. The +x axis is aligned with the streamwise flow. 
The +y axis is horizontal across the tunnel directed towards the left when looking upstream. The +z axis 
is vertical in the tunnel, directed down (see Figure la). 


Table 1. MPF-F111 Parachute parameters and tunnel conditions for all data 


Parachute Parameters 

Tunnel Conditions 

D 0 

1.5816 to (5.189 ft) 

M 

0.106 

So 

1.9554 to 2 (21.048 ft 2 ) 

Q 

736 N/m 2 (15.4 psf) 

R 

3.447 to (11.309 ft) 

Re 

3.7 x 10 6 ** 

m 

0.206 kg 



m a 

0.294 kg 




**Re reference length is D a 


Measurement Apparatus 

The video camera recorded at 29.97 frames per second. Deinterlacing the video signal provided a data frame 
rate of 59.94 frames per second. As this type of data reduction was not intended at the time of the original 
testing, the exact position of the camera was not noted. It is known that the camera was located at the 
same y coordinate in the tunnel as the attachment point of the parachutes. It is also certain that the digital 
still camera (mounted directly beneath the video camera) was mounted at the same z coordinate as the 
parachute confluence point. This would place the video camera approximately 0.152 m (6 inches) above the 
anchor point. A calibration procedure was developed that independently verified this 0.152 m offset. The 
camera’s distance downstream of the parachute is fixed at 5.69 m (224 inches) downstream of the face of the 
capsule backshell. 


Parachute Attitude Extraction 

Figure 2 shows a diagram of the geometry of the parachute in the free-to-rotate test configuration in relation 
to the downstream camera position. The geometry of this setup was used to calculate the pitch and yaw 
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(a) Drag test configuration 



(b) Static aerodynamics test configuration 


Figure 1. TDT Parachute testing setup for drag and static aerodynamics testing. Turntable is on left wall of 
test section as seen looking upstream. 


angles, relative to the flow direction, in the x-z and y-z planes, 9 and ip respectively. These angles were 
measured from each frame of the video footage. The instantaneous orientation of the parachute was defined 
by the position of a line passing through the parachute swivel point and the center of the apex ring. For 
such a non-rigid body, orientation is difficult to define exactly. However, during the static testing by Cruz, 
the parachute was held at angle-of-attack at the apex ring, so these two sets of data should be comparable. 
The apex ring is a metal plug at the apex of the parachute and has a mass of 13.7 grams. This ring was an 
easily identifiable target in the video footage. 

Referring to the geometry in Figure 2, the pitch angle, 9 can be determined from the following equation: 

Rsin9 -j- he ha T he h b 1 1 . . 

l s — Rcos9 l a lb hK 

Here, R is the length of the parachute from the swivel to the apex, h c is the vertical height of the camera 
above the swivel, and l s is the axial distance from the swivel to the camera. The quantity, R, is the axial 
distance from the backshell face to the camera. The constant, K, scales the pixels per meter at the axial 
location of the face of the backshell, and P z is the number of pixels in the video image from the image center 
to the parachute apex in the vertical direction. Equation 1 does not permit an explicit solution for 9, so an 
iterative solution was calculated for each data point. The same relation was calculated for ip except h c = 0 
for the calculations in the x-y plane. Figure 3 shows a captured video image and how the pixels are measured 
for input into equation 1. For each video frame, the distances from the origin, P y and P z , were measured in 
pixels and converted to the pitch and yaw angles, 9 and ip. 
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Figure 2. Geometry of video-parachute setup for orientation extraction from free-to-rotate test 



Figure 3. Video image capture example 


The angles-of-attack and angles-of-sideslip, a and /3, are then determined from the instantaneous pitch and 
yaw angles plus components due to the motion of the parachute as it pivots about its anchor point. 


a = 9 + 


R cp cos 9 L 
V, x 


( 2 ) 


P = ip + 


R cp cos ip 
Vco 


Ip 


( 3 ) 


Here, R cp is the distance from the swivel point to the parachute center-of-pressure (assumed to be the location 
of the parachute center of mass) and T^o is the freestream velocity. For a large freestream velocity, the 8 
and ip terms are small and the total angle-of-attack can be directly measured from the video data. However, 
analysis of this test data indicated that 9 and ip must be included to obtain the total angle-of-attack that 
the parachute canopy experiences when it is free to rotate in the tunnel. 

For the data reduction described later, these angles were combined into a total total angle-of-attack, ax- 
The total angle-of-attack is calculated from the angles, a and fl, as follows: 
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olt = cos 1 (cos a cos (3) 


( 4 ) 


The azimuthal or clock angle, <j>, was also calculated from the pitch and yaw angles to later account for the 
influence of gravity in the equations of motion. The set of total angle-of-attack data and azimuthal angle 
formed the raw data from which the stability coefficients were extracted. 


Calibration 

To determine the pixels-per-meter constant, K , a calibration was performed using video of one of Cruz’ static 
aerodynamic a— sweeps. In Cruz’ test, the full strut mount, holding a parachute, was swept through a range 
of angles-of-attack (—18 to +18 deg). This provided a good set of images to calibrate the pixels-per-meter 
constant as there were several reference points that moved through a known path that was very similar to 
the motions through which the parachute traveled in the free-to- rotate test. Note that for the static test the 
angle-of-attack, a , and the pitch angle in the x-z plane, 9, are identical. 

The geometry of the static aerodynamic test setup is shown in Figure 4. Three points were selected for 
calibration of their observed positions in the video frames to their known positions in the tunnel. These are 
noted in Figure 4. From the geometry of the test setup, the height of the camera in the tunnel, h c , and the 
pixel-per-meter constant, K , were determined from the following three equations which describe the position 
of the three reference points in the video image. 


K 

Ri sin 9 + h c 

Pi 

( 5 ) 

h 

lb — {Ip + R\ cos 6) 

l b K 

K 

f? 2 sin 9 + h c 

P 2 

( 6 ) 

h 

lb — {Ip — R 2 cos 9) 

l b I< 

K 

f ?3 sin 9 — h c 

P 3 

( 7 ) 

h 

lb - {l P + R 3 cos 9) 

lbK 



Figure 4. Static test geometry used for calibration 

Values for h c and K were determined that best fit the set of these equations at each angle-of-attack measured 
in the static test. Using three calibration points added confidence that optical distortion of the image was not 
significant. The values for h c and K were calculated to be 0.1524 m (6.0 in) and 133.86 pix/m (3.40 pix/in) 
respectively. 
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Data Reduction 


Sliding Window Matrix Solution 

To determine the stability coefficients, the time history of the parachute orientation, tracked by the total 
angle-of-attack, ot, and azimuthal angle, </>, was used to extract C m& and C m versus total angle-of-attack. 
For the equations presented here, the derivatives C rn . y and C ma are with respect to total angle-of-attack, 
cut- Other than where needed to account for the influence of gravity, the position in the azimuthal direction 
was not considered in the data reduction. Due to symmetry, the time-averaged values for the static moment 
and damping coefficients orthogonal to olt were assumed to be zero. 

First, the sum of the moments in the olt plane, about the parachute pivot point, is: 


lyy&T PooV^SoDq 


C n 


Do 


-CUT + C ma aT + C mo 


— mg cos 4>R cp cos or = 0 


(8) 


where g is the gravitational acceleration and R cp is the distance from the parachute pivot point to its center 
of pressure (assumed to be at the same point as the center of mass for this analysis). The clock angle, <j>, 
equals zero when t/> is zero and the 0 is positive. The moment of inertia, I yy , about the swivel point is 
dominated by the combined parachute mass and apparent mass (here assumed to be the air within the fully 
inflated canopy) acting as a point mass at a distance, R cp , from the parachute swivel point. The moment of 
inertia about the parachute center of mass is small compared to the (m + m a )R^ p term. Note that the mass, 
m, in the gravity term in Equation 8 is the parachute mass only and does not include the apparent mass. 

In Equation 8, the static aerodynamic moment is represented by the linearized expression: 


M = ^ Poc V^S 0 D 0 [C ma a T + C „ J = \ Poo VlS 0 D 0 C m (9) 

where C ma is the local slope of the static moment coefficient curve with respect to total angle-of-attack and 
Cm 0 is the local intercept of the linearized static moment coefficient curve. Permitting C ma and C mo to 
vary with the local olt allowed for nonlinear moment curves to be extracted from the video data in this data 
reduction. 

For n data points from the video frames, olt was measured and second order finite difference values for olt 
and olt (four-point stencils) were calculated to obtain n-4 versions of equation 8. These were put in matrix 
form to solve for C m6t and C m : 


[ot\ 



r 1 -Pq 

'-'rn.a, 2 Vqo 


C„ 


CL 


OLT 

OLt 

l 


mgRcp r , , 

— [cos 0 cos ar\ 

*yy 


Matrix manipulation yields the expression for the moment coefficients: 


( 10 ) 





T 

- 

r i 


r 

T" 


( 2 \ 

OLT 



OLT 


CtT 


Cm a C mo 

= T /2 c r> ( ; nW mgRcp [cos (j) cos a T \) 

\ Poo r oo OoX'o J 

OLT 



OLT 


OLT 



l 



l 


1 



( 11 ) 


If all data points are used in equation 11, the average values for the stability coefficients across the entire 
data set are obtained. Provided enough data points are available, it is possible to determine the moment 
coefficients as a function of cut- For this analysis, the data were placed in bins of a specified size according 
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to total angle-of- attack. Solving equation 11 for the data in one of these bins finds the moment coefficients 
at a total angle-of-attack somewhere within the bin limits. For this analysis, the coefficients obtained in a 
bin are for the mean total angle-of-attack of all the data in that bin. 

Continuing with this idea, a sliding bin approach was used to obtain the trends of the moment coefficients 
across the total angles-of-attack observed in the video data. A bin size was selected (1.5 deg is typical) and 
placed at the low end of the total angle-of-attack range. For example a bin might contain all points between 
7 and 8.5 deg. Equation 11 was solved for all the data points in that angle-of-attack range. This yields a 
set of stability coefficients at the mean angle-of-attack of the points in the bin. The bin was then slid up 
by an increment, typically smaller than the bin width. For example, the next bin might collect data points 
between 7.1 and 8.6 deg. This new bin contains most of the same points as in the previous bin, losing those 
at angles below 7.1 deg , but adding points between 8.5 and 8.6 deg. Equation 11 was solved again to obtain 
the coefficients at the new mean angle-of-attack and the bin slides over another small increment and so on. 
Proceeding in this manner, values of C m& and C m vs. total angle-of-attack were obtained. 


Validation 

To validate this sliding bin parameter identification, a 2-D angle-of-attack history was generated with known 
moment coefficients. The sliding window data reduction was performed on the angle-of-attack history to 
extract those known coefficients. The angle-of-attack history was generated with a finite difference simulation 
of Equation 8. The coefficients in Equation 8 could be selected arbitrarily, but for this case were based on 
tunnel conditions from the TDT test, the geometry of the MPF F-lll test parachute and the static pitching 
moment data from the static test program (see Table 1). The static measurements of the F-lll parachute 
indicated an unstable trim point at zero degrees angle-of-attack and a stable trim point at olt near 12 
degrees. The static pitching moment data from Cruz’ wind tunnel testing provided a realistic, nonlinear 
pitching moment coefficient curve to test the parameter identification routine. Figure 5 shows this curve 
which was used in the validation simulation. 



Figure 5. Static pitching moment curve used in validation simulation 


To simulate a limit-cycle behavior for the validation simulation only, a simple pitch damping model was 
created where C m& had values of +0.1 from 11 < ar < 13 deg and —0.1 elsewhere. Figure 6 shows the 
validation simulation results and indicates that this damping model with the static pitching moment data 
does produce oscillatory motion that settles into a limit-cycle after an initial displacement from the trim 
angle-of-attack. 

The parameter identification code was then tested by extracting the moment coefficients back out of from 
the simulation results shown in Figure 6. Figure 7a compares the extracted C m with the C m curve used 
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Figure 6. Angle-of-attack time history from validation simulation 


in the simulation for several bin sizes. Figure 7b shows a similar comparison for C md . Figure 7a shows 
the parameter identification methodology is very robust to bin size for extracting the pitching moment. 
The binning technique does have trouble identifying the discontinuity in the simulated pitch damping curve 
although the errors are reduced with smaller bin sizes. A discontinuous pitch damping model is an extreme 
case for the parameter identification technique, but illustrates its limitations. For the parameter identification 
of the parachute video data, a 1.5 deg window was chosen as a compromise between the number of data 
points reduced in each bin and the smoothing introduced by a wide bin. 



(a) Pitching moment (b) Pitch damping 


Figure 7. Aerodynamics extracted from validation simulation against actual simulation values 
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Uncertainty Analysis 


To determine estimates of the uncertainties for the dynamic and static moment coefficient results, the data 
reduction was repeated with perturbed values of the test setup parameters (test apparatus dimensions, 
pixel-to-meter constant, parachute dimensions etc). The data reduction process was repeated 1000 times 
with perturbations placed on the key data reduction parameters. Uncertainty values were specified for each 
parameter and the perturbations were generated randomly with normal distributions within these maximum 
perturbation values. The results of each of the 1000 perturbed cases were collected and averaged. At each 
angle-of-attack the mean and standard deviation of the C m and values were calculated. The 3-standard- 
deviation uncertainty bounds are presented in the results section. The maximum perturbation values are 
listed in Table 2. 


Table 2. Perturbation values for uncertainty analysis 


Parameter 

Uncertainty 

Parameter 

Uncertainty 

Py 

±10 pixels 

I< 

±1.969 pix/m (0.05 pix/in ) 

Pz 

±10 pixels 

m a 

±0.1 kg 

Ibl la,) he 

±0.0254 to (1.0 in) 

R 

±0.1524 to (6.0 in) 


III. Results 


Parachute Attitude Measurements 

Approximately 79 seconds of video data were analyzed, yielding 4761 data points from which pitching moment 
and pitch damping information was extracted. Figure 8 shows a segment (for clarity) of the parachute’s 
orientation history captured from the video data. The data points are the raw positions extracted from the 
video. The curve is a local average of the nearest three points along the time history. Comparing the raw 
data with observations of the video footage, it appeared that the raw data captures some of the motion of the 
parachute disk relative to the band. This may be exacerbated by the presence of the apex ring. Recall that 
the apex ring is a significant amount of mass at the apex of the parachute (6.6% of the entire parachute/riser 
mass). The video footage suggests that the relative motion between the disk and band adds high frequency 
fluctuations to the bulk motion of the parachute. As the apex ring was used in the video analysis to identify 
the parachute position, this relative motion was captured in the raw data. In an attempt to better represent 
the parachute motion due to aerodynamic forces, three-point averaged data was used for the parameter 
identification. Note that the strut extending from the wall to support the backshell model and swivel mount 
created a wake, disrupting the flow downstream toward the parachute. The parachute tended not to travel 
through the wake region of the strut. This is apparent in Figure 8 where 6 « 0 deg and ip < 0 deg. Any 
disruption to the aerodynamic behavior of the parachute caused by the strut’s wake was not quantified in 
the data reduction. 


Parameter Identification 

Pitching Moment (olt plane) 

The 1.5 deg sliding window data reduction results for pitching moment are shown in Figure 9. These data 
show the parachute to be stable and trims at a total angle-of-attack of 10 deg. The uncertainty bounds 
determined from 1000 perturbed cases, shows an uncertainty on the trim angle of approximately ±1.25 deg. 
The parameters that dominate the uncertainties are the parachute length, pixel-per-meter constant and 
the apparent mass of air in the parachute. Perturbations of the parachute length and the pixel-per-meter 
constant tend to shift the trim point, while perturbing the apparent mass (which affects the moment of 
inertia) changes the slope of the pitching moment curve only. 
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Figure 9. Pitching moment coefficient from 1.5 deg sliding window parameter identification compared to static 
tunnel data by Cruz. 4 Moment reference point is at the parachute anchor point. 

Figure 9 also shows the static test results (with 3-standard-deviation confidence limits) from Cruz. 4 Note 
that the data determined from the video images show the parachute trims 2 degrees shallower than does the 
static test data. This appears to be a real difference between the parachute’s apex-unconstrained behavior 
and behavior while forced to an angle-of-attack. The difference in trim point is greater than what can be 
accounted for by uncertainties in each test. This disagreement is not unexpected as significant deformation of 
the parachute canopy was observed at high angles-of-attack during the static testing. The pitching moment 
slopes of both sets of data agree fairly well at large angles-of-attack. At total angles-of-attack below the 
trim point the video-derived pitch damping curve becomes more shallow again agreeing more closely with 
the slope of the static test data. 
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Pitch Damping (ax frame) 


Figure 10 shows the pitch damping results for 1.5 and 3.0 deg windows. For such chaotic motion of the 
parachute, it had been hoped that a trend be found for pitch damping as a function of total angle-of-attack. 
The deformation of the parachute canopy and shedding vortices produce forces and moments that cause 
dramatic changes in the direction and speed of the parachute, seemingly at random. Even with 79 seconds of 
data, the 1.5 deg window data reduction does not produce a smooth C mA curve. To add data points to each 
bin (at the expense of angle-of-attack resolution) a 3.0 deg window was used. The 3.0 deg window results 
do fit within the smaller window data and vary more continuously. Note that a window which encompasses 
the entire data set would yield one value for C m<i across this entire total angle-of-attack range. Assuming 
the parachute has reached a quasi-steady state (e.g., a limit cycle), this would require C m& to be zero for 
the full data set. 



Figure 10. Pitch damping coefficient from parameter identification using 1.5 and 3.0 deg sliding windows. 
Moment reference point is at the parachute anchor point. 


The uncertainty analysis helps understand the apparent trend in pitch damping. Along with the uncertainty 
bounds, the mean C m& values of the 1000 perturbed cases are plotted in Figure 10. This mean curve bisects 
the 1.5 and 3.0 deg window data and shows a trend with angle-of-attack. At the trim angle-of-attack, the 
parachute is dynamically unstable (C m<i >0). At low angles, the parachute becomes dynamically stable. 
Above 11.5 deg and below 7 deg total angle-of-attack, the amount of data points drops off significantly and 
the uncertainty bounds diverge and become quite large. Sets of approximately 2000 equations of the form 
shown in Equation 8 were used to solve for the stability coefficients (see Equation 11) near the trim angle- 
of-attack, while sets of approximately 100 equations were used at the high and low angle limits of the data. 
The variation in available data points at each angle-of-attack accounts for the variation in uncertainties. 
The dynamic stability variation across the entire total angle-of-attack range are small compared to the 
uncertainties. Therefore, definitive comments about the pitch damping of the parachute are not possible 
from these data. However, dynamic instability at the trim angle-of-attack coupled with regions of dynamic 
stability away from the static trim angle is a plausible mechanism for the bounded total angle-of-attack 
range observed in the free-to-rotate video. Dynamic instability at the trim angle would tend to drive the 
parachute from equilibrium, while dynamically stable regions away from the trim angle would tend to damp 
out large motions and the statically stable parachute would return to the trim point. 

The bin-averaged pitch damping results are insufficient to replicate the complete behavior of the parachute. 
The instantaneous pitch damping values can deviate significantly from the bin-averaged data points shown 
in Figure 10. A more complete model of the parachute behavior may require random components be added 
to the pitch and yaw damping and static stability coefficients as a function of time. Also, the uncertainties 
introduced by the selection of bin size remain a concern, although there is fairly good agreement between 
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the 1.5 and 3.0 degree window data. Furthermore, the static and dynamic data extracted here provide for 
no mechanism that would introduce azimuthal motions. A model capable of predicting azimuthal motions 
that vary in direction and magnitude, while settling to a quasi steady state is required to fully describe the 
parachute motion. Such a model can not be extracted by the methods developed here, although it should 
be possible to expand on the principles used in this effort to do so. 

IV. Conclusions 

A method for extracting pitching moment versus total angle-of-attack from video data of a parachute that 
is free to rotate about its anchor point has been developed. Validation work using a simulated oscillating 
parachute first showed that this parameter identification technique can extract nonlinear static and dynamic 
pitching moment information from the motion of a 2-D simulation. Next, these parameter identification 
techniques extracted a pitching moment curve from 79 seconds of video data. These data show the parachute 
to be statically stable at a total angle-of-attack, in agreement with static wind tunnel tests. The trim angle- 
of-attack, however, was identified to be two degrees shallower than was measured in the static tunnel tests. 
A study of bin size showed that the 79 seconds of video data were sufficient to capture “bin-converged” static 
pitching moment characteristics of the parachute. 

Although the erratic motion of the parachute suggests that the pitch damping of this parachute fluctuates 
significantly from instant to instant, the parameter identification results indicate that, averaged over time, 
the pitch damping varies with angle-of-attack. These data suggest that the parachute is dynamically unstable 
in the ar plane at the static trim angle-of-attack and becomes dynamically stable as the parachute’s attitude 
departs from this trim angle. Bin-converged solutions were not possible for the dynamic derivatives although 
the mean curve from uncertainty perturbation cases fits the different bin size results well. Uncertainty 
analysis indicates that the details of the observed trends are small compared to the uncertainty bounds. 

The coefficients extracted with this data reduction technique offer additional information to complement 
the earlier static testing and free-to-rotate testing. However, these data are still insufficient to fully model 
the observed motion of the MPF-F111 disk-gap-band parachute. There is no functionality in the static or 
dynamic coefficients to produce changes in the direction of motion in the azimuthal direction or non-periodic 
motion in the otr plane. Variations from the time-averaged values of the static pitch and yaw moments and 
pitch and yaw damping (due to vortex shedding, canopy shape fluctuations, etc.) may all be significant, 
but cannot be extracted by the methods presented here. Future work should expand the application of 
this new technique to attempt to quantify the frequency and time variation of the stability coefficients from 
the time averaged curves identified here. Extracting pitch and yaw stability coefficients with the complete 
equations of motion may be possible by improving upon this technique, but would require much more data. 
The technique as presented took advantage of an azimuthal symmetry assumption to extract the stability 
coefficients as a function of olt only. 

This work was done as a proof-of-concept activity. Future parachute tests, with more forethought, may 
use these techniques to extract static and dynamic stability data with more accuracy. Well placed and 
clearly identifiable markers located on the parachute, careful calibration of the camera system and accurate 
measurement of the test setup geometry should yield more accurate and precise results. 
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